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We discuss the phenomenon of long-distance entanglement in the ground state of quantum spin 
models, its use in high-fidelity and robust quantum communication, and its realization in many- 
body systems of ultracold atoms in optical lattices and in arrays of coupled optical cavities. We 
investigate XX quantum spin models on one-dimensional lattices with open ends and different 
patterns of site-dependent interaction couplings, singling out two general settings: Patterns that 
allow for perfect long-distance entanglement (LDE) in the ground state of the system, namely such 
that the end-to-end entanglement remains finite in the thermodynamic limit, and patterns of quasi 
long-distance entanglement (QLDE) in the ground state of the system, namely, such such that the 
end-to-end entanglement vanishes with a very slow power-law decay as the length of the spin chain is 
increased. We discuss physical realizations of these models in ensembles of ultracold bosonic atoms 
loaded in optical lattices. We show how, using either suitably engineered super-lattice structures 
or exploiting the presence of edge impurities in lattices with single periodicity, it is possible to 
realize models endowed with nonvanishing LDE or QLDE. We then study how to realize models 
that optimize the robustness of QLDE at finite temperature and in the presence of imperfections 
using suitably engineered arrays of coupled optical cavities. For both cases the numerical estimates 
of the end-to-end entanglement in the actual physical systems are thoroughly compared with the 
analytical results obtained for the spin model systems. We finally introduce LDE-based schemes of 
long-distance quantum teleportation in linear arrays of coupled cavities and show that they allow 
for high-fidelity and high success rates even at moderately high temperatures. 



PACS numbers: 03.67.Hk, 03.67.Mn, 75.10.Pq 
I. INTRODUCTION 

Quantum entanglement plays a crucial role in many ar- 
eas of quantum information science [l| , including, among 
others, quantum cryptography and secure quantum key 
distribution and quantum communication [jj. For 
this reason, much work has been dedicated to single out 
ways to produce useful entanglement for efficient and ro- 
bust implementation of quantum information tasks. The 
most natural way to create entanglement between two 
or more constituents of a quantum information device 
appears to be by means of direct interactions because, 
intuitively, large amounts of entanglement should be as- 
sociated to the presence of strong quantum correlations. 
However, from a general quantum informatic perspective, 
an even more desirable goal is to envisage ways to pro- 
duce large amounts of entanglement shared between dis- 
tant and generally not directly interacting constituents. 
Along the way to accomplish this task, a first step has 
been realized by introducing the concept of localizable 
entanglement [J], namely the rate of entanglement that 
can be concentrated on a pair of arbitrarily distant con- 
stituents by performing optimal local measurements onto 
the remainder of the system. More recently, it has been 
shown that the ground state of some spin models with fi- 
nite correlation length defined on one-dimensional chains 
with open ends can support large values of the end-to-end 
entanglement between the initial and final points of the 



chain, insensitive or only very weakly sensitive to the size 
of the system^, 0, 0, HI- This type of entanglement has 
thus been dubbed Long-Distance Entanglement (LDE). 
In this approach, the guiding principle is to look for sys- 
tems whose GS can support intrinsically large amounts 
of (long-distance) entanglement, even at moderately high 
temperatures and in the presence of bulk imperfections, 
without the need for performing operations and measure- 
ments or for continuous dynamical controls and adjust- 
ments of the couplings. Clearly, the property of LDE 
would be very appealing in the quest to design quantum 
information devices, that perform tasks efficiently and 
at the same time are sufficiently robust against decoher- 
ence, based on (properly engineered) many-body systems 
of condensed matter that could possibly be realized with 
currently available technologies, or technologies that may 
be in view in the near future. 

The scope of the present work is to investigate var- 
ious theoretical aspects of the phenomenon of LDE in 
quantum many-body physics, to introduce possible con- 
crete experimental implementations with currently avail- 
able many-body systems of atomic physics and quantum 
optics, and to discuss simple schemes of LDE-based ef- 
ficient and robust quantum communication. Concerning 
the theoretical modeling, here and in the following we 
will consider as paradigmatic test-beds some very simple 
isotropic quantum spin models on open chains with XX- 
type interactions of site-dependent strength. For suit- 
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able choices of the sets of site-dependent couplings, these 
models can show large values of the entanglement shared 
by the end spins of the chain 0,11]. Notwithstanding their 
apparent simplicity, they are a very powerful and useful 
tool in order to investigate, identify, classify, and summa- 
rize the different possible types of LDE, and the prop- 
erties associated to different patterns of site-dependent 
interaction strengths. Moreover, these different classes 
of XX quantum spin models with site-dependent inter- 
action couplings turn out to be very well suited to assess 
the experimentally relevant problem of the vanishing of 
the energy gap as a function of the size of the system, 
and to identify, unambiguously, the optimal range of pa- 
rameters compatible with the largest achievable value of 
LDE at fixed values of the temperature and of the length 
of the chain. 

After establishing the general theoretical framework, 
we will analyze some possible realizations based on two 
many-body systems of current experimental interest and 
that are moreover evolving and being developed at a fast 
pace, namely systems of ultracold neutral atoms in opti- 
cal lattices Q and arrays of coupled cavities [1JJ] . These 
two systems share some common features. They can be 
seen as ensembles of local structures that interact with 
each other via the exchange of bosonic particles, atoms 
in the case of optical lattices, and photons in the case of 
coupled cavities. Both systems allow to simulate quan- 
tum spin models by effectively reducing the dimension 
of the local state space, and by realizing spin-spin inter- 
actions via the inter-site hopping amplitudes. However, 
while the hopping amplitude in an array of coupled opti- 
cal microcavities can be tuned at will, at least in princi- 
ple, either by adjusting the inter-cavity overlaps and/or 
by tuning the cavities' fundamental physical parameters 
(Rabi couplings, quality factors, etc.) of every single cav- 
ity (single-site addressing), this task appears to be much 
more challenging in optical lattices, for which the inter- 
site properties are indissolubly woven with the frequency 
of the external laser potentials, and hence site-dependent 
interactions between the effective spins can be engineered 
either by introducing atomic impurities and/or by design- 
ing super-lattice structures. We will discuss merits and 
disadvantages of both types of systems in the quest for 
the experimental demonstration of LDE and its use in the 
realization of efficient, robust, and high-fidelity quantum 
information tasks. 

The paper is organized as follows: In Section II we 
introduce the general XX spin chain Hamiltonian with 
arbitrary site-dependent couplings and review the meth- 
ods to determine its spectrum, its ground-state proper- 
ties, and the end-to-end concurrence. In Section III we 
investigate the properties of LDE for different patterns 
of site-dependent interaction couplings, singling out two 
different general situations: Patterns that allow for per- 
fect LDE, namely such that the end-to-end entanglement 
remains finite in the limit of diverging size of the sys- 
tem, and quasi-LDE (or imperfect LDE), i.e. such that 
the end-to-end entanglement vanishes with a very slow 



power-law decay as the length of the chain increases. 
In Section IV we discuss how to realize or, better, how 
to simulate the models analyzed in Section III with en- 
sembles of ultracold bosonic atoms loaded in optical lat- 
tices. In the literature there exist several proposals to 
implement LDE-free spin models, either with uniform 
[Til H~3. f]~3j] or random [3] spin-spin couplings, using op- 
tical lattices in various settings. Here we show how, using 
either suitably engineered super-lattice structures or ex- 
ploiting the presence of edge impurities in lattices with 
single periodicity, it is possible to realize models endowed 
with nonvanishing LDE in the ground state. In section 
V we investigate how to realize spin models supporting 
LDE by the use of suitable array structures of coupled 
cavities. In fact, these systems are being intensively stud- 
ied in relation to their ability to realize/simulate collec- 
tive phenomena typical of strongly correlated systems of 
condensed matter [H EE EE LLl EE El- Concerning 
the simulation of spin models, we show how, at least in 
principle, this type of system allows to simulate straight- 
forwardly XX models with arbitrary patterns of site- 
dependent couplings and different types of LDE. In both 
sections numerical estimates of the end-to-end entangle- 
ment in the actual physical systems are thoroughly com- 
pared with the results obtained in the simulated spin 
model systems. In Section VI we discuss the implemen- 
tation of an LDE-based scheme of long-distance quan- 
tum teleportation exploiting open linear arrays of cou- 
pled cavities. In particular, we show that this scheme 
allows for high-fidelity teleportation even at moderately 
high temperatures and in the presence of noise. Finally, 
in the conclusions we summarize our findings and discuss 
some outlooks on possible future developments along this 
line of research. 



II. GENERAL STRUCTURE OF XX QUANTUM 
SPIN MODELS ON OPEN CHAINS 

As anticipated in the introduction, in this Section we 
discuss the properties of XX quantum spin-1/2 models 
defined on one-dimensional lattices with open ends, with 
specific patterns of nearest-neighbor interactions. Such 
models are all special instances of the general XX Hamil- 
tonian 

N-l 

H XX = J2 Ji (SfSf +1 + SfSf +1 ) , (1) 
»=i 

where Ji is the interaction strength between spins at 
nearest neighboring sites i and i + 1, Sf is the spin-1/2 
operator defined at site i, and N is the total number of 
sites (spins) or, equivalently, the length of the chain. The 
spectrum of this Hamiltonian can be determined exactly 
by a straightforward generalization of the methods first 
discussed by Lieb, Schultz, and Mattis [20]. The first step 
in the procedure is to perform a Jordan- Wigner transfer- 
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where Sj 1 — Sj ± iSj . As a result, the Hamiltonian ([T]) 
is mapped in the free fermion Hamiltonian 
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fermionic creation (annihilation) operators, one for each 
site of the lattice, and the adjacency matrix M reads 
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Given the GS of the system, we are interested in the 
evaluation of the end-to-end entanglement in the reduced 
two-qubit state between the first spin at site 1 and the 
last spin at site N. One thus needs to determine the 
spin-spin concurrence (entanglement of formation) be- 
tween the end points of the chain. This quantity can 
be computed exactly for any two-qubit state (pure or 
mixed), thanks to the celebrated Wootters formula [22| . 
and the task is left to obtain its explicit expression in 
the reduced state of the two end-point spins. To this 
purpose, we need to calculate explicitly all the possible 
forms of two-point correlations in the GS. The Hamilto- 
nian (fTJ) is symmetric under rotations of the spins around 
the z-axis, so that the only nonvanishing correlations are 
{SfSp = (SfSV), (SfSp and (Sf). In the absence of 
external magnetic fields, n- rotations around the x and 
y axes are symmetries of the model, which additionally 
implies (Sf) = at every site. Thanks to the aforemen- 
tioned symmetries, the two-point reduced density matrix 
Pi j , obtained by tracing the ground-state density matrix 
of the entire system over all spins except the ones at sites 
{i,j}, has the form 

Pi tj = - + (SfSJ) (a x ®<j x +a y ® a y ) + (S*S*)a z <g> a\ 

(5) 

where <r x ' v < z are the Pauli matrices and (•) is the GS 
average at temperature T = 0, or the thermal one at fi- 
nite temperature (5 = (^T)" 1 with respect to the Gibbs 
state p = e~ l3H Z~ 1 . We are interested in the case in 
which i = 1 and j = N are the two end points of the 
chain. In this instance, we have 
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where Af = Yl!i=i c \ c i IS *he total number operator. 
Applying Wick's theorem and taking into account that 
1/2, we obtain 
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Setting x = (c\cn) we see that the end-to-end reduced 
density matrix depends uniquely on this parameter. For 
reduced states of the form ([5]) the end to end concurrence 
C\.n is readily computed [231 ] . In our case, we have 
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The concurrence is thus nonvanishing for |x| > (y/2 — 
l)/2 ~ 0.207, and it reaches the maximum value C\ t N — 
1 for \x\ -►1/2. 

It is natural to expect that the existence of a strong 
quantum correlation between the two end spins of the 
chain can be conveniently exploited for performing tasks 
in quantum information, in particular considering tele- 
portation schemes. In the standard quantum teleporta- 
tion protocol, two parties A and B share a maximally 
entangled state (Bell state). Party A holds also a third 
qubit, whose unknown state is to be teleported. If the 
two end points of our XX chain share a highly entangled 
state, that in some limit may even be asymptotically close 
to a Bell state, they can be identified as the two parties, 
sender and receiver, for a long-distance, high-fidelity tele- 
portation protocol. The efficiency of a quantum channel 
in teleporting an unknown state is quantified by the fi- 
delity / between the output and the input states, aver- 
aged over all input realizations. The fidelity depends on 
the actual properties of the entangled resource p\ : L (Cfr. 
Eq. shared by the end spins of the chain. In fact, it 
has been demonstrated that the optimal fidelity depends 
only on the "fully entangled fraction" Ff u u , according to 
the formula / = (2F fuU + l)/3 p3|. The fully entangled 
fraction is defined as the fidelity between the resource 
pi t L and a maximally entangled state, maximized over 
all possible maximally entangled states. For states of 
the form Eq. (0 it can be easily computed, and reads 
2 [IH. The associated teleportation 
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and, taking into account Eq. {8J, from Eq. ([9]) we obtain 
the expression of the concurrence C\ t N as a function of 
the fidelity: 



C hN = 2max<j 0, -/ - 1 



(10) 



This expression highlights the crucial interplay between 
entanglement and efficiency in quantum information pro- 
tocols. In fact, due to the high symmetry of states of 
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the form Eq. ([5]), a nonvanishing entanglement implies 
a nonclassical teleportation fidelity exceeding the classi- 
cal threshold 2/3, and viceversa. Moreover the relation 
between the maximum fidelity that can be obtained in 
a quantum teleportation between the end points of the 
spin chain and the end-to-end Eq. () 1 0(1 provides a sim- 
ple route to probe, via experimental implementations of 
quantum teleportation protocols, the actual values of the 
generated long-distance entanglement C\^m- 

All the physical information about model ([3]) can now 
be obtained by diagonalizing the one-body matrix M. 
Let £fc be the eigenvector of M with eigenvalue where 
k is a quasi-momentum label. Passing to new fermionic 
operators via the transformation Cj = Efc£fc °k, the 
Hamiltonian takes the form 



H X x = ^2 Afccjcfc 



(11) 



The evaluation of the two-point correlation x is finally 
straightforward: 
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where x depends on the set of the couplings J; as well as 
on the temperature T, whenever T =/= 0. 

In conclusion, the full evaluation of the end-to-end en- 
tanglement in XX quantum spin models defined on open 
chains, with arbitrary patterns of site-dependent nearest- 
neighbor interactions, is strictly associated to the com- 
plete diagonalization of the M matrix in Eq.([3|). While 
in some particular cases it is possible to find the exact 
analytic expressions of the eigenvalues and of the eigen- 
vectors of M @, Hi| , in the most general case the problem 
has to be solved by exact numerical diagonalization. 



III. END-TO-END ENTANGLEMENT 
PROPERTIES OF XX SPIN CHAINS 

Having reviewed the basic definitions and mathemat- 
ical tools needed to analyze the end-to-end entangle- 
ment properties in XX spin chains with arbitrary site- 
dependent interaction strengths, we need to look for 
those sets of couplings that allow for large values of LDE, 
possibly robust against thermal decoherence and bulk im- 
perfections, the required entanglement properties. The 
basic idea for concentrating a large amount of entangle- 
ment on the two end spins of the chain is inspired by 
the observation that for some frustration-free systems, 
like the antiferromagnetic Heisenberg model H @ , and 
some XX-type models 0, [1], the total energy is mini- 
mized by states tending to form a global singlet. This 
fact, together with t he p roperty of monogamy for bi- 
partite entanglement [27| . is the primary cause for the 



phenomenon of LDE. Namely, if one selects a set of in- 
teraction couplings in the Hamiltonian Eq.|T|) that are 
such to forbid the end spins to entangle with any other 
constituents, then, by monogamy, strong quantum corre- 
lations will develop between them. 



A. Perfect Long Distance Entanglement (LDE) 

In analogy with the case of the dimcrized antiferromag- 
netic Heisenberg chain, that was the system in which the 
phenomenon of LDE was first discovered [B[ , we begin by 
reviewing the properties of the XX spin chain with bonds 
of alternating strengths, i.e. a model in which each weak 
bond (Jfc-i = AJ, with A << 1) between neighboring 
spins at sites is followed by a much stronger one ( Jk = J) 
between the successive pair of neighboring spins. This 
model has been introduced and investigated in Ref. 0]. 
If the chain is formed of an even number of sites and 
the two end points are coupled to their nearest neigh- 
bors with a weak bond, this model fulfills exactly the 
requirements needed in order to establish perfect LDE 
between the end spins of the chain. Namely, any spin 
in the chain but the end ones is subject to a weak (left) 
interaction and to a strong (right) interaction. Hence, 
every spin in the bulk tends to be maximally entangled 
with its strongly interacting neighbor. Each of the two 
end spins is subject only to a weak bond interaction, and 
is therefore excluded from the arrangement in pairs of 
nearest-neighbor singlets due to the monogamy property 
of bipartite entanglement [27| ■ Furthermore, because the 
GS of the model tends to be a collection of singlets, the 
two end spin spins are effectively coupled and develop 
a nonvanishing LDE mediated by the remainder of the 
chain, notwithstanding the fact that they do not interact 
directly in any way. 

The above qualitative reasoning can be made precise 
and quantitative. The Hamiltonian for the XX spin 
model with alternating weak and strong interactions is 



H = 'k £ K 1 + A ) + t- 1 ^ 1 - A )i ( Wi + WO ■ 

»=1 

Following the approach described in the previous section, 
one can determine the end-to-end concurrence for differ- 
ent values of A and different lengths of the chain 0], as 
shown in Fig. [TJ We can see that as the weak coupling A 
decrease below a threshold value A c = 0.765 [3], a nonva- 
nishing concurrence begins to develop. For a fixed length 
of the chain, decreasing A still further (and thus increas- 
ing the difference between the alternating couplings) al- 
lows the creation of dimers between pairs of strongly in- 
teracting neighboring spins. The dimers, in turn, are 
very weakly coupled to each other. Hence, because the 
two end spins of the chain interact very weakly with their 
nearest neighbors, dimerization and monogamy of entan- 
glement force the creation of a strong quantum correla- 
tion between the end points, in analogy with what hap- 
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Figure 1: End-to-end concurrence Ci,n for the 1-D XX spin 
model with alternating couplings described by the Hamilto- 
nian Eq |13l plotted as a function of the number of lattice 
sites N (length of the chain), for different value of the weak 
coupling A. From bottom to top: magenta empty squares: 
A = 0.7; brown empty circles: A = 0.6; orange inverted tri- 
angles: A = 0.5; blue triangles: A = 0.4; green diamonds: 
A = 0.3; red full squares: A = 0.2; black full circles: A = 0.1. 




Figure 3: Energy Gap AE(J) between the ground and the 
first excited state, in units of the coupling energy J, for the 
1-D XX spin model with alternating couplings described by 
the Hamiltonian Eq. (|13|l , plotted as a function of the length 
N of the chain for different values of the weak coupling A. 
From top to bottom: magenta empty squares: A = 0.7; brown 
empty circles: A = 0.6; orange inverted triangles: A = 0.5; 
blue triangles: A = 0.4; green diamonds: A = 0.3; red full 
squares: A = 0.2; black full circles: A = 0.1. 
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Figure 2: End-to-end entanglement as a function of the length 
of the chain for the model with noisy couplings described by 
Eq. (I15|l , where x% is a random variable uniformly distributed 
in the interval {— x,x}> an d X = 0.2. The violet area cor- 
responding to A = 0.2, and the blue one corresponding to 
A = 0.4, are the domains in which the different random sam- 
ples may fall. They have been obtained with 10 4 independent 
samples for each value N of the length of the chain. 



pens in the Heisenberg case [5|. On the other hand, at 
fixed A, the end-to-end concurrence Ci.jv a s a function 
of the size N of the chain converges very rapidly to its 
thermodynamic-limit saturation value: 

Ci i00 -2max|o,i-A 2 + ^| . (14) 

Next, it is important to verify that the phenomenon of 
strong end to end entanglement survives in the presence 
of disorder. Indeed, in Fig. [5] we represent the behavior 
of the end-to-end entanglement for different sample re- 
alizations of the following Hamiltonian with disordered 



couplings: 

iV-l 

H = 2 £ [(l+A)+(-ir'(l-A)](l+Xi) (ST + S?S? +1 ) , 
t=l 

(15) 

where \i is a random variable uniformly distributed in 
the interval {—A 7 , A 7 }- As we may note, the presence of 
noise affects bi-directionally the entanglement properties 
of the system, and the end to end concurrence is non- 
vanishing even in the presence of a noise with a rel- 
ative weight of 20%. For nearly half of the samples, 
noise increases the end to end concurrence. Moreover, 
Fig. [3 shows that the stronger the end-to-end entangle- 
ment without noise, the weaker the relative effect of the 
noise. 

It would then seem that in order to obtain a phys- 
ical system with strong LDE it should be sufficient to 
engineer some concrete device that in appropriate limits 
simulates/realizes an XX spin chain chain with alternat- 
ing weak and strong couplings. Unfortunately, this is not 
the case in realistic situations. As shown in Fig. [3l the 
phenomenon of perfect LDE is strictly associated to an 
energy gap between ground and first excited states that 
vanishes exponentially fast as a function of the length 
of the chain. Therefore, as soon as the system is at 
temperatures comparable with the gap AE, the equilib- 
rium Gibbs state results in an incoherent superposition 
of singlet and triplet states in which the LDE vanishes 
0, 0] ■ Consequently, in order to realize a nonvanishing 
entanglement between the two end points, even in short 
chains, it would be necessary to reach temperatures fan- 
tastically close to absolute zero, something that is well 
beyond reach in any foreseeable future, and in any case 
of no use for practical applications. 

One might still consider realizing the request of per- 
fect LDE by looking at modifications or alternatives to 
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Figure 4: Schematic drawing of XX spin chains endowed with 
perfect LDE in the GS. Upper drawing: Model Eq.([T3J, with 
alternating weak (A < 1) and strong couplings along the entire 
chain. Lower drawing: Model Eq. (|16p , with left and right end 
regions of alternating weak and strong couplings and a central 
(bulk) region of uniform strong couplings. 



the pattern of perfectly alternating weak and strong cou- 
plings. Such a pattern has been introduced with the 
aim of suppressing the formation of sizeable entangle- 
ment between the end points and the remaining spins 
of the lattice. However, it is quite clear that this goal 
can be achieved considering different spatial interaction 
patterns that have the effect of isolating the bulk of the 
chain from the end points. This observation motivates 
the introduction of the following class of Hamiltonians: 
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(16) 



Models Eq. (fT51) and Eq. (jTHJ) with their spatial patterns of 
site-dependent couplings are sketched pictorially in Fig. 
|H The choice of the spatial pattern of the couplings 
that characterizes model Eq. pH)) allows to preserve per- 
fect LDE between the end points of the XX spin chain, 
at the same time drastically increasing the energy gap 
between the GS and the first excited states. In practice, 
it amounts to replace in the bulk of the chain, i.e. be- 
tween a reference N-th spin and the N — N-th one, the 
alternating pattern of weak and strong couplings with a 
uniform nearest-neighbor interaction. The results for the 
end-to-end entanglement as a function of the length of 
the chain, for different values of A and at a fixed value of 
the ratio N/N are reported in Fig. O In the thermody- 
namic limit, the model described by Eq. (|16p exhibits a 
LDE behaviour analogous to that of the model with fully 
alternating weak and strong couplings described by Eq. 




Figure 5: End-to-end concurrence Ci,n for an XX spin chain 
with alternating couplings and a central region with uniform 
interactions, described by the Hamiltonian Eq. (|16|l . plotted as 
a function of the length N of the chain, for different values of 
the weak coupling A and a fixed value of the ratio N/N = 1/4. 
From bottom to top: magenta empty squares: A = 0.7; brown 
empty circles: A = 0.6; orange inverted triangles: A = 0.5; 
blue triangles: A = 0.4; green diamonds: A = 0.3; red full 
squares: A = 0.2; black full circles: A = 0.1. 




Figure 6: Energy gap, in units of the coupling energy J, be- 
tween the ground and the first exited state for an XX spin 
chain with alternating couplings and a central region with 
uniform interactions, described by the Hamiltonian Eq, (|16[l . 
The gap AE(J) is plotted as a function of the length TV of 
the chain, for different values of the weak coupling A and at a 
fixed value of the ratio N/N — 1/4. From top to bottom: ma- 
genta empty squares: A = 0.7; brown empty circles: A = 0.6; 
orange inverted triangles: A = 0.5; blue triangles: A = 0.4; 
green diamonds: A = 0.3; red full squares: A = 0.2; black full 
circles: A = 0.1. 



(|13[) . when the size of the chain is increased at a fixed, 
constant value of the ratio N/N. The only difference, as 
can be seen by comparing Fig. Q] with Fig. is that 
in the second case the end-to-end concurrence reaches 
the saturation value more slowly, with a speed that de- 
creases when decreasing the ratio N/N. While the es- 
sential features of the LDE behaviour are shared also by 
the modified model with uniform central interactions, the 
behaviour of the energy gap between the ground and the 
first excited states is instead strongly affected, as shown 
in Fig. [5] As we can see from Fig. [5J the energy gap again 
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decreases as the length of the chain is increased, as in the 
previous case, and it attains values that forbid possible 
experimental implementations of LDE at finite temper- 
ature (although higher than the corresponding ones for 
the model with perfectly alternating couplings). Com- 
paring Fig. [6] with Fig. [3] and taking into account that 
in the present case N — N/4 the behavior of the energy 
gap does not depend on the total length N of the chain, 
but rather only on the length 2N = N/2 of that part of 
the chain characterized by alternating pairings. This de- 
pendence of the energy gap on the length of that part of 
the chain endowed with alternate couplings is confirmed 
by the analysis of systems defined at different values of 
N. This fact suggests to investigate further possibilities. 
Indeed, one may think of bringing the modification im- 
plemented in Eq. (|16p to the extreme, in order to realize 
non negligible values of LDE at small, but experimentally 
feasible temperatures. 

B. Quasi Long Distance Entanglement (QLDE) 

As we have seen at the end of the previous subsection, 
reducing the portion of an XX chain with alternating 
patterns of interaction does not reduce the entanglement 
shared by the end points and increases by several orders 
of magnitude the energy gap between the ground and the 
first excited states, yet, unfortunately, it does not allow 
LDE to survive except at unrealistically low tempera- 
tures. Therefore, in order to realize a nonvanishing LDE 
at low, but realistically attainable temperatures, we con- 
sider taking the limit of Eq. (|16[) to a model with uniform 
nearest-neighbor interactions for all pairs of spins but for 
the two end-points, that are connected to the rest of the 
chain with a weak bond. Indeed, spin systems allowing 
for strong end-to-end correlations should be character- 
ized by interactions between the end points and their 
nearest neighbors that are always smaller than the inter- 
actions in the bulk of the chain. Otherwise, if the system 
does not meet this criterion, the end points would be- 
come strongly entangled with their neighbors, excluding, 
due to the monogamy constraints the possibility of 
LDE. Hence, we consider an open XX spin chain formed 
by N — 2 spins with uniform coupling strengths, plus two 
weakly interacting probes placed at the two end points. 
Such a model is described by the Hamiltonian 

N-2 

H = [SfSi+i + (17) 

_L T \ ( Q x Q x _L Qy Qy _L Q x ox , qV qV \ 
-\-JA {Di 2 T &x &2 ' N-l°N °N-l D N) ' 

where < A < 1. The evolution of the end-to-end con- 
currence shared by the end points as a function of the 
length of the chain for different values of the weak cou- 
pling A is reported in FigJT] At variance with what oc- 
curred in the models with patterns of alternating cou- 
plings studied in the previous subsection, in the case 
of models with uniform bulk interactions and weak end 




Figure 7: End-to-end concurrence Ci,n for an XX spin chain 
with uniform bulk couplings and two weak end probes de- 
scribed by the Hamiltonian Eq. (|17)) . plotted as a function of 
the length N of the chain for different values of the weak 
coupling A. From bottom to top: orange inverted triangles: 
A = 0.2; blue triangles: A = 0.16; green diamonds: A = 0.12; 
red full squares: A = 0.08; black full circles: A = 0.04. 



probes the GS end-to-end entanglement is sensitive to 
the size of the chain. It decreases slowly as the length 
of the chain grows and vanishes at a critical value of the 
length, which depends on the value of the weak coupling 
A. This behavior of the GS end-to-end entanglement is 
due to the fact that the region responsible of preventing 
the formation of quantum correlations between the end 
points and the bulk of the chain is fixed at two lattice 
spacings, and does not grow as the length of the chain is 
increased. However, by choosing sufficiently small values 
of the weak coupling A, it is always possible to realize 
a GS with nonvanishing LDE in chains of arbitrary fi- 
nite length. Hence it is natural to name this behaviour 
as quasi LDE (QLDE), a phenomenon that differs from 
perfect LDE, as the latter occurs only when the system 
admits a non vanishing GS end-to-end entanglement that 
attains its maximum (saturation) value in the thermody- 
namic limit. 

In addition to the QLDE nature of the end-to-end en- 
tanglement in the GS, models with weak end bonds pos- 
sess an energy gap between the GS and the first excited 
states that exhibits a very different behavior compared 
to the case of models characterized by LDE. As one can 
see from Fig. [FJ the energy gap does not exhibit an 
exponential decay with the size of the chain as in the 
cases of Fig. [3] and Fig. [6l This fact, together with 
the possibility to choose appropriately the value of A to 
ensure a non vanishing QLDE shared by the end points 
in chains of arbitrary finite length, opens the possibil- 
ity for the experimental realization of QLDE at low but 
realistically achievable temperatures with concrete phys- 
ical systems engineered in configurations suitable for the 
realization/simulation of model Eq. (fl~7|) . 

Even if the models introduced in the previous subsec- 
tion and the one described by Eq. (fT7| show different 
behaviors of the end to end entanglement, they share the 
property that it is robust in the presence of noise. Indeed, 
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Figure 8: Energy Gap AE(J), in units of the coupling energy 
J, between the GS and the first excited states for an XX spin 
chain with weak end probes, described by the Hamiltonian 
Eq. (|17p . plotted as a function of the length N of the chain for 
different values of the weak coupling A. From top to bottom: 
orange inverted triangles: A = 0.2; blue triangles: A = 0.16; 
green diamonds: A = 0.12; red full squares: A = 0.08; black 
full circles: A = 0.04. 
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Figure 9: End-to-end entanglement as a function of the 
length of the chain for the model with disordered couplings 
described by Eq. (|18[) where Xi is a random variable uniformly 
distributed in the interval {— X, x} with X = 0.2. The violet 
area corresponding to A = 0.2, and the blue one correspond- 
ing to A = 0.4, are the domains in which the different random 
samples may fall. They have been obtained with 10 4 inde- 
pendent samples for each value TV of the length of the chain. 



Fig. [9] displays the end-to-end entanglement for different 
realizations of the model with disordered couplings de- 
scribed by the following Hamiltonian: 



N-2 



H = J^Xi [S*S? +1 + SfSf +1 \ 

z=2 

J(X + XN-i)(S x N _ 1 S% + S y N _ 



(18) 



where \i l& a random variable uniformly distributed in 
the interval {— x,x}- As with models endowed with per- 
fect LDE, also in the case of QLDE the presence of noise 
affects the end to end concurrence bi-directionally, and 
the entanglement remains nonvanishing even when the 
relative weight of the noise is up to 20%. Also in the 



Figure 10: Schematic drawing of XX spin chains endowed 
with QLDE in the GS. Upper drawing: Model Eq.(TJ7J), with 
weak end bonds A < 1 and uniform unit couplings in the 
bulk. Lower drawing: A-^i model Eq. (|19jl . with weak end 
bonds A < 1, strong near-end bonds n > 1, and uniform unit 
couplings in the bulk. 



case of models associated to QLDE, about half of the 
random sample realizations happen to increase the value 
of the concurrence between the end points. Moreover, 
the disruptive effect of noise is strongly reduced for large 
values of the end-to-end entanglement. 

The role played by the two weak interactions in the 
model Eq. (TT7|) is to force the second (N— 1-th) spin of the 
chain to get heavily entangled with the third (N — 2-th) 
spin, thus forcing the first and last spins of the chain, due 
to entanglement monogamy and the fact that this con- 
figuration is energetically favorable, to develop a nonvan- 
ishing quantum correlation in the GS. This effect can be 
enhanced further by increasing the interaction between 
the second and the third spin, and correspondingly be- 
tween the N — 1-th and the N — 2-th, well above the 
reference value between neighboring spins in the bulk of 
the chain. In this case the model Hamiltonian reads: 



iV-3 



i=3 



(19) 



qV qV \ qx qx 

+v (sf sf + sf sf + s 



qy 

N- 



QX 

2 D N- 



qy qy 

J N — 2 J N- 



Models Eq. (fTT|) and Eq. p^|) with their spatial patterns of 
site-dependent couplings are sketched pictorially in Fig. 
[Tm In Eq.fUni) we have A < 1 and [i > 1, so that AJ 
(weak end bond) < J (uniform bulk interaction) < [iJ 
(strong near-end bond). The increased efficiency granted 
by model Eq. (|19|) for the process of creating a nonvanish- 
ing QLDE between the two end points of the chain can 
be appreciated in Fig. [TT1 where the model (fTF)) with 
simple weak end bonds A = 0.1, /i = 1, is compared 
with model (fT9"|) with the same value of A and different 
values of the strong near end bonds \i > 1. In the fol- 
lowing, Hamiltonians (| T T[) realizing simple QLDE in the 
GS will be referred to as A models, while Hamiltonians 
(fT9|) that realize enhanced QLDE in the GS will be re- 
ferred to as X-fi models. It is not surprising that the 
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Figure 11: End-to-end concurrence Ci,jv for the X-/J, spin 
chain described by Hamiltonian Eq. (|19[l . plotted as a function 
of the length N of the chain, for a fixed value of the weak 
coupling A = 0.1, and different values of the strong coupling 
fi. From top down: orange inverted triangles: fi — 5; blue 
triangles: /i = 4; green diamonds: /i = 3; red full squares 
fi = 2. The lowest-lying line (black full circles) is the one 
with fi — 1, corresponding to the A model Eq. (|17p . and is 
drawn for comparison. 
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Figure 13: Two-dimensional contour plot of the end-to-end 
concurrence Ci,jv as a function of the weak and strong cou- 
plings A and fj,, for a 1-D A-/i spin model (|19[) on an open linear 
chain of N = 10 spins at a reduced temperature T/J = 0.05. 
The color map relative to the values of the concurrence is 
reported on the left. Values vary from (violet) to 1 (red). 
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Figure 12: Energy Gap AE(J) (in logarithmic scale), in units 
of the uniform bulk coupling J, between the GS and the first 
exited states for the A-/i spin chain described by Hamiltonian 
Eq. (|19[l . plotted as a function of the length N of the chain, for 
a fixed value of the weak coupling A = 0.1, and different values 
of the strong coupling fj,. From bottom up: orange inverted 
triangles: /i = 5; blue triangles: fi — 4; green diamonds: 
fj, = 3; red full squares fj, = 2. The line with [i — 1 (black full 
circles) corresponds to the A model Eq. (|17[) . and is drawn for 
comparison. 



strong enhancement brought by the X-/J, model to the 
QLDE that can be accommodated in the GS is obtained 
at the cost of a (relatively moderate) trade-off with the 
behaviour of the energy gap. In Fig. [T^] we have com- 
pared the behaviour of the energy gap, as a function of 
the size of the system, for the A model and various X—fi 
models with the same value of A and different values of 
fi. At this point, the question naturally arises of optimiz- 
ing the parameters of the Hamiltonian in order to single 
out the maximum possible value of the QLDE compat- 
ible with an energy gap sufficiently large to warrant a 
concrete physical realizability at finite temperature. Ob- 
viously, the process of optimization of the Hamiltonian 



parameters does not provide a unique general answer. 
Rather, the results will vary, depending on the length of 
the chain that one wants/needs to consider as well as on 
the minimum working temperature that is fixed by the 
external experimental conditions. In the following we will 
always consider the models of enhanced QLDE described 
by Hamiltonian §T§\i , that for [i = 1 reduce to the models 
(fTT|) of simple QLDE. A typical result of the optimization 
process for a X—fi spin chain of finite length is shown in 
Fig. [T3]in the case N = 10. It is important to notice that 
too high values of fi and/or too low values of A, contrary 
to naive intuition do not help, because they lead cither 
to an effective separation of the end points from the rest 
of the chain and/or to a much too small energy gap com- 
pared to the energy amplitude of thermal excitations. In 
both cases, the associated QLDE vanishes. Moreover, the 
fact that very large values of the end-to-end entanglement 
are always associated with small energy gaps implies, as 
the figure shows, that the optimization process is never 
trivial (especially in determining the regions of matching 
values for the couplings), and therefore it represents a 
key point in the discussion of realistic physical systems 
able to simulate/realize models endowed with QLDE, a 
subject that we investigate in the following sections. 



As with the other cases considered previously, also the 
A — model is robust against noise affecting the coupling 
amplitudes. This resilience is shown in Fig. 1141 display- 
ing the end-to-end entanglement for different realizations 
of the model with disordered couplings described by the 
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Figure 14: End-to-end entanglement as a function of the 
length of the chain for the model with disordered couplings 
described by Eq. p0[) where Xi is a random variable uniformly 
distributed in the interval {— x, X.} with X = 0.2. The violet 
area, corresponding to A = 0.1 and [i = 2.0 and the blue 
area, corresponding to A = 0.1 and fj, = 5.0, are the domains 
in which the different random samples may fall. They have 
been obtained with 10 4 independent samples for each value 
of the length N of the chain. 

Hamiltonian: 

N-3 

H = J ^2 Xi [S?S* +1 + SfSf +1 ] (20) 

t=3 

+X X i (S?Sf + SfSl) 

+ MX2 (s%si + sist) 

Here, as in the previous cases, Xi 1S a random variable 
uniformly distributed in the interval {— 

IV. END-TO-END ENTANGLEMENT IN 
OPTICAL LATTICES 

Until now we have analyzed in some detail the ques- 
tion of identifying specific instances of open XX spin 
chains with different patterns of site-dependent interac- 
tions that possess GSs endowed with large values of the 
entanglement shared by the end points and an energy gap 
compatible with physical realizations at finite tempera- 
ture. Here and in the following we investigate concrete 
instances of atomic and quantum optical systems that 
may allow the experimental demonstration of LDE and 
QLDE. 

We first consider systems of ultracold neutral atoms 
loaded on a one-dimensional optical lattice generated by 
interfering laser beams [9(. The dynamics of such sys- 
tems is usually well described, depending on the nature 
of the atoms present in the ensemble, by a Bose-Hubbard 
model 28] or by a Fermionic Hubbard model, or by mix- 
tures of atoms of different species, such as Bose-Bose (29j . 
Fermi- Fermi [3(j , and Bose-Fermi mixtures [3l[ . Regard- 
less of the specific cases realized, the parameters in the 



various Hamiltonians can be easily controlled by tuning 
both the intensity and the frequency of the laser beams. 
This fact, i.e. that all the Hamiltonian parameters can 
be manipulated and modified with a high degree of con- 
trol, makes optical lattices of particular interest for the 
simulation/realization of models of interacting quantum 
systems. The mapping of Hubbard-type models into 
XX spin models, in appropriate regimes of the parame- 
ters, was originally discussed in the fundamental work of 
Fisher et al. [32j where the Bose-Hubbard model was 
introduced. In recent years, several works have been 
dedicated to the simulation of spin models using Hub- 
bard Hamiltonians realized in optical lattices. However, 
in most of these works the spin-spin interactions are usu- 
ally taken to be either uniform [13, or completely 
random [lil ]. 

Let us consider a 1-D optical chain loaded with single- 
species ultracold bosonic atoms. The system is described 
by the 1-D Bose-Hubbard Hamiltonian [28[ that can be 
written as 

11 N N ^ 
H=-J2rn(ni-l)-tJ2(Hb i+ i + h.c.) , (21) 

i=l i=l 

where b\ (bi) is the creation (annihilation) operator of a 
bosonic atom at the i-th site of the chain, m = b\bi is the 
number operator at site i, and N is the number of sites in 
the chain. The parameter U denotes the strength of the 
local on-site repulsion, while t is the hopping amplitude 
between adjacent sites. These Hamiltonian parameters 
depend on the quantities characterizing the external pe- 
riodic optical field, taking into account only the lowest 
vibrational states for every minimum of the periodic po- 
tential [28j [. 

The simplest approach to simulate XX spin models 
starting with the bosonic Hubbard Hamiltonian Eq. (l2~T1) 
relies on the fact that the local Fock states with either 
one or no particle per site differ from all other local Fock 
states in that they are the only two local levels whose 
energy is not dependent on the repulsion. Therefore, in 
the strong interaction regime (U ^> t), they are automat- 
ically separated from the others, and hence they can play 
the role local spin-1/2 states. In the framework of this 
approximation, the hopping term is naturally mapped in 
an interaction of the XX type with an amplitude equal 
to two times the hopping amplitude. However, in an opti- 
cal lattice, the hopping amplitude between adjacent sites, 
depends both on the frequency and the intensity of the 
lasers that are quantities common to each site. Hence, 
the interaction couplings between adjacent spins in the 
ensuing XX model are site- independent. 

However, site-dependent spin-spin couplings can be 
easily engineered by introducing local fields on properly 
selected sites of the lattice. Consider a local field on a 
given site, say the fc-th one, whose amplitude is larger 
than the hopping amplitudes but weaker than the on- 
site repulsion U. In this situation occupation of the on- 
site single-particle Fock state is energetically unfavorable. 
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Therefore, if an atom in k + 1-th or in k — 1-th site hops 
onto the A;-th site, it is immediately pushed away. This 
mechanism realizes either a self-interaction at site k or an 
effective interaction between the k + 1-th and the k — 1-th 
site. This simple reasoning can be put in a quantitative 
form resorting to degenerate perturbation theory in pow- 
ers of the hopping amplitude. Let k — 1, k, and k + 1 
be three adjacent sites, let Ek ^ the local field at site 
k, and such that t <C <C U. In second-order per- 
turbation, discarding all the excited states with energy 
proportional to U and , each pair term in Eq. (|21[) maps 
into the following spin-spin interaction: 




It is important to note that in these spin-spin interac- 
tion terms the dependence on the fc-th site is removed. 
Hence, the introduction of site-dependent spin-spin inter- 
action strengths is associated to a reduction of the sites of 
the optical lattice. Eq. ([22|) is the basic ingredient in the 
realization of XX spin chains with site-dependent inter- 
actions and nonvanishing end-to-end entanglement. The 
successive task is to determine suitably engineered op- 
tical lattices with the appropriate local field dynamics. 
For instance, one possible way to realize the XX spin 
model with alternating weak and strong couplings Eg J13I 
is to introduce an optical super-lattice potential obtained 
through the interference of two sets of laser beams of dif- 
ferent amplitude strengths. The set of stronger beams 
realizes the lattice structure while the set of weaker ones 
realizes and modulates the presence of a local field at 
each lattice site. If, for instance, the set of weak beams 
is tuned at a wavelength three times as large as that of 
the set of strong ones, one realizes a system in which ev- 
ery site over three is characterized by a strong local field, 
while the remaining two sites experience a practically 
vanishing local field, as schematically illustrate in Fig. 
1151 Resorting to degenerate second-order perturbation 
theory in powers of t, and taking into account Eq. (f!?2)) . 
the optical lattices Hamiltonian is mapped in the spin 
Hamiltonian 

i— 1 ^ ' 

where e > t is the maximum amplitude of the local field, 
and N' = 2 A/3 is the number of spins in the chain that, 
as already noticed, does not coincide with the number of 
sites in the physical optical lattice chain. The last sum 
of terms represents an overall energy offset that can al- 
ways be reabsorbed in the definition of the zero of the 
energy and has no dynamical effect. Fixing the total 




Figure 15: Scheme of the optical super-lattice potential ob- 
tained by combining two appropriate sets of stationary laser 
beams. It realizes an effective XX spin chain with alternat- 
ing weak and strong interactions that allows for LDE between 
the two end points. The first set of laser beams generates an 
optical potential (black line) that realize an optical lattice. 
On the contrary, the weaker optical potential of the second 
sets, represented by the red line, creates a local field at each 
site of the optical lattice. The local field is strongly enhanced 
on every site over three. 



magnetization at 0, corresponds to fixing the number of 
atoms in the optical lattice at n = N'/2 = N/3. In this 
case we obtain that in the strong interaction regime the 
optical super-lattice simulates the 1-D XX spin chain 
with alternating weak and strong couplings Eq |131 ex- 
actly with A = t/e. It is important to verify the sound- 
ness of the theoretical mapping, obtained in second-order 
perturbation theory, by direct exact numerical compar- 
ison between the end-to-end entanglement properties of 
the spin model Eq[l3] and of the original bosonic system. 
In Fig. [T|)]we report the end-to-end entanglement in the 
GS of a system of ultracold bosonic atoms loaded on the 
above-described optical super-lattice, with N = 12 sites 
and n = N/3 = 4 atoms, corresponding, according to 
the our previous analysis, to a 1-D XX spin chain with 
N' = 2N/3 = 8 sites with alternating weak and strong 
couplings. As measures of entanglement, we consider the 
concurrence [22| and the logarithmic negativity 33] . The 
concurrence quantifies the GS entanglement between the 
two end points seen as effective spins, obtained by trac- 
ing out all Fock states but the empty state and the state 
with single-atom occupation. The logarithmic negativity 
is used to measure the actual GS entanglement of the full 
bosonic system between the two end sites of the optical 
super-lattice. Fig. [16] shows that for small values of the 
on-site repulsion U there is no end-to-end entanglement 
regardless of the value of the local field e. Viceversa, in 
the strong interaction regime the system develops a size- 
able end-to-end entanglement that always approximates 
very closely the end-to-end entanglement of the XX spin 
chain with alternating weak and strong couplings, prov- 
ing that appropriate optical lattice or super-lattice struc- 
tures can simulate/realize efficiently XX spin chains with 
LDE entanglement. An exception is constituted by the 
limiting case e = t. In this situation the strong local 
field condition e > t is violated, A = 1 in the correspond- 
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Figure 16: End-to-end entanglement in the GS of a super- 
lattice system of TV = 12 sites and n = 4 atoms, with strong 
local field e on 2 — nd, 5 — t/i, 8 — th, and 11 — th site, plotted 
as a function of the on-site repulsion U in units of the hopping 
amplitude t. The dashed lines represent the logarithmic nega- 
tivity, measuring the end-to-end entanglement in the original 
Bose-Hubbard atomic system, while the solid lines stand for 
the concurrence, that measures the end-to-end entanglement 
in the effective XX spin model with alternating weak and 
strong couplings. The different curves correspond to different 
values of the strong local field e. From bottom up: orange: 
e = t; green: e = 3t; blue: e = 5t; red: e = It; black: e = lOi. 

ing spin Hamiltonian, and thus the end-to-end spin-spin 
concurrence vanishes, while the end-to-end logarithmic 
negativity in the GS of the original optical super-lattice 
system is definitely nonvanishing. 

We now turn to the question of simulating spin models 
with QLDE rather than pure LDE. One problem in sim- 
ulating the Hamiltonians described in Subsection IIII Bl 
with atomic systems in optical potential structures is 
that QLDE models do not have a periodic pattern of 
site-dependent interactions and thus are difficult to ap- 
proximate with super-lattices, as done for LDE models. 
Spin models with QLDE structures can be simulated by 
addressing the degrees of freedom of the individual lat- 
tice sites with pre-selection techniques that have recently 
been discussed for different purposes [H| , [H| . If these 
methods, or similar ones, could be implemented success- 
fully in realistic experimental conditions, it would be pos- 
sible to realize optical-lattice structures such that only 
few pre-selected sites are affected by a nonvanishing local 
field. Consider the following Bose-Hubbard Hamiltonian 

H = -^2ni(ni-l)-t^2(blb i+ i + h.c.) 

i=l i=l 

+e(n 2 + njv-i) . (24) 

The presence of a local field e at the second and N — 1-th 
sites allows to simulate, when t < e < U, the XX spin 
model with weak end probes described by Eq. (fTT|) . In 
Fig. (fl~7|) wc report the GS end-to-end entanglement for 
the Hamiltonian Eq. (|2~4"|) defined on an optical lattice 
of N — 10 sites and loaded with n — 4 bosonic atoms, 
and compare it with the GS end-to-end entanglement of 
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Figure 17: End-to-end entanglement in the GS of the Bose- 
Hubbard Hamiltonian Eq. (|24p . defined on an optical lattice 
of TV = 10 sites with strong local field on the 2 — nd and on 
the 9 — th site, plotted as a function of the on-site repulsion U 
in units of the hopping amplitude t. The dashed lines corre- 
spond to the logarithmic negativity measuring the end-to-end 
entanglement in the original Bose-Hubbard atomic system, 
while the solid lines stand for the concurrence, that measures 
the end-to-end entanglement in the effective XX spin model 
with weak end probes. The different curves correspond to 
different values of the local field e. From bottom up: green: 
e — 3t; blue: e = 6i; red: e = 9i; black: e — 12t. 
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Figure 18: Energy gap between the GS and the first excited 
states of the Bose-Hubbard Hamiltonian Eq. (I24|l , defined on 
an optical lattice of TV = 10 sites with strong local field on 
the 2 — nd and on the 9 — th site, plotted as a function of 
the on-site repulsion U in units of the hopping amplitude t. 
The different curves correspond to different values of the local 
field e. From bottom up: green: e = 3t; blue: s = 6t; red: 
e = 9i; black: e = 12t. 

an XX model with weak end probes Eq. (fTT|) defined on 
a chain with 8 spins. We see that the former is an ex- 
cellent approximation to the latter in the limit of large 
on-site repulsion U. We have fixed the total number of 
bosonic atoms to (N-2)/2 to satisfy the half-filling con- 
dition. This corresponds again to impose a vanishing 
magnetization in the effective spin chain. In the strong- 
coupling regime, the value of the end-to-end entangle- 
ment grows as the local field e is increased, as should 
be expected from the results of the discussion in Subsec- 
tion HTlBl taking into account Eq. (|22p. However, as 
discussed in the previous section, the significance of sim- 
ulating/realizing spin models with QLDE using, for in- 
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stance, optical lattice Hamiltonian structures as the ones 
described by Eq. is that the energy gap can be suf- 
ficiently large to ensure that a nonvanishing end-to-end 
entanglement is still retained at very low but realistically 
reachable working temperatures. The two requirements 
of large values of QLDE and a large energy gap are partly 
conflicting. The behavior of the energy gap is reported 
in Fig[THl Comparing Figs [T71 and IT51 we see that the op- 
timal compromise, that allows to obtain both a sizeable 
QLDE and a sizeable energy gap, takes place in initial 
side of the strong-coupling regime, at moderately large 
values of the on-site repulsion U. 

V. END-TO-END ENTANGLEMENT IN 
ARRAYS OF COUPLED OPTICAL CAVITIES 

Recently, hybrid atom-optical systems of coupled cav- 
ity arrays (CCAs) have been intensively studied in re- 
lation to their ability to simulate collective phenomena 
typical of strongly correlated systems [bsL [lrl fl7l. [Til. fl9l| . 
In the present section we will discuss in detail some recent 
findings about the possibility of exploiting appropri- 
ately engineered arrays of coupled optical cavities to re- 
alize XX open spin chains sustaining LDE or QLDE. Be- 
sides the extremely high controllability and the straight- 
forward addressability of single constituents, arrays of 
coupled cavities also allow in principle a great degree of 
flexibility in their design and geometry [1 01 ] . Therefore, 
they should be tested as natural candidates for the re- 
alization of spatially extended communication networks 
and scalable computation devices. 

Consider a linear CCA with open ends, consisting of 
N cavities. The dynamics of a single constituent of the 
array doped with a single two-level atom is well described 
by the Jaynes-Cummings Hamiltonian [3(| 

H k = uja\a k + uj'\e k }(e k \ + ga\\g k ) (e k \ + g\e k }(a k \a k , 

(25) 

where a k (a k ) is the annihilation (creation) operator of 
photons with energy lo in the fc-th cavity, \g k ) and \e k ) are 
respectively the ground and excited atomic states, sepa- 
rated by the gap a/, and g is the photon-atom coupling 
strength. The local Hamiltonian Eq. (|25|) is immediately 
diagonalized in the basis of dressed photonic and atomic 
excitations (polaritons) : 

|0*> = k)|0 fe ); (26) 
\n+ k ) = cos8 n \g k )\n k ) + sin0 n |e fc )|(n - l) fc ) n > 1 ; 
\n—h) = sm8 n \g k )\n k ) - cos9 n \e k )\(n - l) k ) n>l, 

where 9 n is given by tan 29 n — —g^/n/A and A = to' — lj 
is the atom-light detuning. Each polariton is character- 
ized by an energy equal to 

£o = 0; e„± =nw± ^Jng 2 + A 2 . (27) 

When u = V 'g 2 + A 2 the GS of Eq. $Z5§ becomes two- 
fold degenerate, resulting in a superposition of |0fc) and 
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Figure 19: Representation of the energy levels for a cavity 
with lj — \/ g 2 + A 2 . The GS is two-fold degenerate, and the 
energy gap £2- prevents the occupancy of the higher energy 
levels, thus realizing an effective two-level system. 

|1— k), see Fig. [19] If both the atom-cavity interaction 
energy and the working temperature are small compared 
to £2- = 2y ' g 2 + A 2 — \/2g 2 + A 2 , one may neglect all 
the local polaritonic states but \% k ) and |1— This sit- 
uation defines a local two-level system. Adjacent cavities 
can be easily coupled either by photon hopping or via 
wave guides of different dielectric and conducting prop- 
erties. The wave function overlap of two adjacent cavities 
introduces the associated tunneling elements, so that the 
total Hamiltonian of the CCA reads: 

N N-l 

H cca = ^Hk- ^ Jk(a[a k+1 + a) k+1 a k ) . (28) 

k k 

Each hopping amplitude J k depends strongly on both 
the geometry of the cavities and the actual overlap be- 
tween adjacent cavities. If the maximum value among 
all the couplings {J k } is much smaller than the energy 
of the first excited state: max{Jfc} <C £2-, then the to- 
tal Hamiltonian Eq. (j2"5)) can be mapped in a spin-1/2 
model of the XX type with site-dependent couplings of 
the form Eq. (fl~9l) . where the state |0fe) (|1— k)) plays the 
role of I Ik) (| Tfc))- The mapping to an open-end A-/i 
linear spin chain sustaining QLDE in the GS is then re- 
alized, e.g., by simply tuning the distance between the 
end- and next-to-end sites of the array, as showed in 
Fig. [20] The A-/i XX spin chain thus can be realized 
starting from an equispaced CCA with site-independent 
nearest- neighbor bulk coupling amplitude Jf,, and then 
engineering appropriately the positions of the 2 — nd and 
of the (N — 1) — th cavities, that are placed closer to 
their neighbors in the bulk and farther away from the 
end points of the array. This shift lowers below unity the 
reduced coupling between the end points of the array and 
their next-to-end neighbors (Ji/Jb — Jjv-i/ 'Jb = A < 1), 
and increases above unity the reduced coupling between 
the next-to-end sites and their neighbors in the bulk 
{J21 1 Jb = JN-2/Jb = Li > 1), effectively realizing the 
A-^i model of QLDE. Obviously, adjusting the hopping 
rate by spacing the cavities closer to the end points is 
not the only way to realize the A — /x model in an array 
of optical cavities. Indeed, it can be technically rather 
challenging. An alternative way to tune the inter-cavity 
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Figure 20: Scheme of an array of coupled optical cavities re- 
alizing a X-fj, linear spin chain described in Eq. (|19p . Dark 
green: Area covered by the wave functions associated to each 
site of the array. The two next-to-end sites (light green cir- 
cles) are symmetrically displaced with respect to a situation 
of perfectly equispaced sites, and drawn closer to their neigh- 
bors in the bulk (black circles). As a consequence, the overlap 
(black area) between the wave functions of these two cavities 
and their neighbors in the bulk is larger than the would-be 
reference (unit) overlap in an equispaced array. At the same 
time, the overlap between the end sites of the array (red and 
blue circles) and the next-to-end sites is reduced proportion- 
ally compared to an equispaced array. 



hopping amplitudes, probably much more feasible with 
currently available technologies, is to dope the cavities 
with few-level atoms and use Raman transitions and the 
atomic Lambda level structures to tune the couplings be- 
tween two cavities in situ, by using external laser drives 
as described, e.g., in Ref. [13 |. 

In the next section we will discuss the properties of 
QLDE, realized using CCAs, in connection with the im- 
plementation of tasks of quantum information science. 
In particular, we will discuss how the A-/z model and the 
associated GS QLDE can be exploited to realize high- 
fidelity, long-distance quantum teleportation protocols. 



VI. APPLICATIONS: QLDE AND 
LONG-DISTANCE QUANTUM 
TELEPORTATION IN CCAS 

We now proceed to illustrate that CCAs in the A- 
fi configuration allow for long-distance and high-fidelity 
quantum communication in realistic conditions and at 
moderately high temperatures. In Fig. [5T] we report 
the fidelity of teleportation F max [24[ as a function of 
the reduced couplings A and \i for different temperatures. 
Remarkably, Fig. [5U shows the existence of a rather high 
critical temperature of teleportation for CCAs realizing 
a A— fx spin chain. The region of the physical parameters 
compatible with a nonclassical fidelity F max > 2/3 is 
progressively reduced with increasing temperature, until 
it disappears at T c rs 0.13J;,. Similar behaviors are ob- 
served for longer CCAs, with T c slowly decreasing with 
the length of the array. For instance, for an array of 
N = 36 cavities in the A— fi configuration, the critical 
temperature of transition to bona fide quantum telepor- 




Figure 21: (Color online) Fidelity of teleportation F max in 
a X-fi configuration, by exact numerical diagonalization, for 
a CCA of N = 12 cavities as a function of the couplings 
A = Ji/Jb and fi — 3^1 Jb at different temperatures T/Jb. 
Left panel: T/J b = 0.005. Right panel: T/J b = 0.01. F max 
varies between 0.5 (violet) and 1 (red). Horizontal dashed 
line: Classical threshold Ff nax = 2/3. 

tation is T c w 0.11J&. 

A formidable obstacle to the concrete realization of 
working quantum teleportation devices is performing the 
projection over a Bell state, in order for the sender to 
teleport a quantum state faithfully to the receiver. In 
fact, in the framework of condensed matter devices there 
hardly exist quantities, easily available in current and 
foreseeable experiments, that admit as eigenstates any 
two-qubit Bell states. To avoid this problem different 
schemes for quantum information transfer have been de- 
veloped, for instance, for what concerns entanglement 
transfer, in Ref. [37|- In the following, we will discuss 
a simple and concrete scheme for long-distance, high- 
fidelity quantum teleportation in A— fJ, CCAs that realizes 
Bell-state projections indirectly, by matching together 
free evolutions and local measurements of easily control- 
lable experimental quantities [H, HH, • We first illus- 
trate it in the simplest case of two cavities at zero temper- 
ature, with the first cavity accessible by the sender and 
the second one by the receiver. The sender has access 
also to a third cavity, the "0" cavity, that is decoupled 
from the rest of the chain, and stores the state to be tele- 
ported \tp) = a\ To}+/3| lo)- The decoupling between the 
— th cavity and the rest of the chain can be achieved 
removing the degeneracy among |0)o and |1— )o and tak- 
ing |eo — Ex— | 3> Jo- The total system is initially in the 
state 

|*(0)> = ±(a\ To)+/3| lo»(| Ti)l ia>+| Ta». (29) 

At t = the state begins to evolve and if Jo 3> J\ one 
has: 

!*(*)} = i[«lTo)|Ti)U2>+j8Uo)Ui)|ta) (30) 

I To) | ii)(acoB(J t)| h)-^|3sm(Jot)\ j 2 » 
| | )| Ti>Hasin(J *)| h) +/3cos(J t)| | 2 ))] . 

If at time t — 7r/(4Jo) Alice measures the local magne- 
tizations (Sq, Sf) in the first two cavities, she will find 
with probability 1/2 that the teleported state is the im- 
age of \tp) under a local rotation. The value 1/2 for the 
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Figure 22: (Color online) Fidelity of teleportation / in a X-fi 
CCA channel of N = 12 cavities, for a generic input \ip) = 
a \ To) + P\ lo), as a function of \a\ at different temperatures. 
Black line: T = 0.001 Jt,; Red line: T = 0.003 J;,; Green line: 
T = 0.004J b ; Blue line: T = 0.005Ji,; Magenta line: T = 
0. 007 J b . Here A = 0.5, ^ = 4.0, and v = 50. Horizontal 
dashed line: Classical threshold f c = 2/3. 

probability stems from the fact that any simultaneous 
eigenstate of Sq and Sf can be obtained with equal prob- 
ability but one may discard the case in which the total 
magnetization is equal to ±1. Realizing a local rotation 
of ±7r/2 around Sf, with the sign depending on the re- 
sult of the measurement that the sender communicates 
classically to the receiver, the latter recovers the original 
state \ip) with unit fidelity. The simple protocol described 
above can be immediately extended to X-fi CCAs of any 
size, at finite temperature, and removing the constraint 
J\ <C Jo- By resorting again to exact diagonalization, in 
Fig. [22] we report the behavior of the average fidelity of 
teleportation /, as a function of |a| of the state |</?), in the 
case of an array of N — 12 cavities, with v = Jo/Jb = 50 
and for different temperatures. Also in the non-ideal case 
the teleportation protocol has probability 1/2 of success. 
The fidelity depends on the input state, with a maxi- 
mum for inputs with |a| = |/3| = 1/V2 and a minimum 
for inputs with |cv| = 0.1. The fidelity remains above 
0.95 for all values of \a\ at moderately low temperatures 
(T~10- 3 J 6 ). 

VII. CONCLUSIONS AND OUTLOOK 

In conclusion, we have introduced some classes of quan- 
tum spin models defined on open 1-D lattices that are 
characterized by non-perturbative GSs with nonvanish- 
ing LDE or QLDE. In particular, end-to-end QLDE is 
found to be strongly resilient to thermal decoherence and 
can be efficiently achieved with a minimal set of local ac- 
tions on the end- and near-end couplings in linear open 
CCAs. We have showed that QLDE-based CCAs allow 
for a simple quasi-deterministic protocol of long-distance, 
high-fidelity quantum teleportation that yields a high 
rate of success without direct Bell measurements and pro- 
jections over Bell states. In detail, we have studied mod- 
els of XX quantum spin chains with nearest-neighbor 
interactions, discussing two types of spatial patterns of 



the coupling strengths. In the case of chains with bonds 
of alternating strengths, we have shown that the exact 
GS possesses a nonvanishing LDE between the end spins 
of the chain, independently of the size of the system and 
asymptotically close to unity (maximal entanglement) in 
the limit of exact dimerization. This system is therefore 
perfectly suited for bona fide long-distance quantum tele- 
portation with ideal fidelity at zero temperature. How- 
ever, the limiting maximal values of the fidelity are ob- 
tained at the cost of introducing an energy gap above 
the GS that vanishes exponentially with the size of the 
system. Therefore, this model is de facto useless for ef- 
ficient quantum teleportation at finite temperature. We 
have then discussed another class of XX open spin chains 
with uniform bulk interactions and small end bonds. In 
this case, we have shown that for sufficiently small values 
of the end couplings, the GS of the system possesses a 
sizeable QLDE between the end spins of the chain, that 
at fixed size of the system is asymptotically close to unity 
(maximal entanglement) in the limit of vanishingly small 
end bonds. However, QLDE decreases, albeit slowly, as 
the size of the system is increased. Models of QLDE 
can be improved further by introducing strong near-end 
bonds that enhance the value of the end-to-end entan- 
glement (A-/i model). An interesting feature of QLDE 
models is that the energy gap above the GS vanishes 
only algebraically, as the first power of the inverse of the 
size of the system. Therefore, in principle, spin chains 
endowed with the QLDE property can be exploited as 
quantum channel for teleportation with nonclassical fi- 
delity at finite temperature. 

A relevant issue concerned the behavior of these sys- 
tems in the presence of disorder. This question is im- 
portant especially in view of possible experimental im- 
plementations in which the couplings can be engineered 
only within a certain accuracy. Since naturally the effect 
of disorder is that of localizing eigenstates, we expected 
LDE to be more robust than QLDE against disorder. 
In fact, in the LDE scenario the states responsible for 
the end-to-end entanglement is already localized at the 
borders whereas localization is only approximate in the 
QLDE case. This conjecture has in fact been confirmed 
by thorough exact numerical samplings, that also show 
that the relative weight of imperfections that can be tol- 
erated by the systems, maintaining large and useful val- 
ues of LDE and QLDE, is rather high, up to 20%, and 
more. Moreover, for about always half of the spatial pat- 
terns and values of the couplings, disorder can even en- 
hance the end-to-end entanglement. It is at the moment 
unclear how this picture extends to finite temperature. 
Another interesting open problem worth further study 
is to assess the existence and the possible location of a 
crossover between perfect LDE and prima facie QLDE. 

Demonstrating experimentally (Q)LDE and (Q)LDE- 
based efficient long-distance teleportation and state 
transfer will be a first crucial preliminary test in order 
to proceed with integrated devices of quantum informa- 
tion science, combining atom-optical systems such as op- 
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tical lattices or systems of trapped ions with solid-state 
based systems, such as those of circuit Cavity Quantum 
Electrodynamics. 
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